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Abstract — We provide a simple and accurate analytical model 
for multi-cell infrastructure IEEE 802.11 WLANs. Our model 
applies if the cell radius, R, is much smaller than the carrier 
sensing range, R cs . We argue that, the condition R cs >> R 
is likely to hold in a dense deployment of Access Points (APs) 
where, for every client or station (STA), there is an AP very 
close to the STA such that the STA can associate with the AP 
at a high physical rate. We develop a scalable cell level model 
for such WLANs with saturated AP and STA queues as well 
as for TCP-controlled long file downloads. The accuracy of our 
model is demonstrated by comparison with ns-2 simulations. We 
also demonstrate how our analytical model could be applied in 
conjunction with a Learning Automata (LA) algorithm for opti- 
mal channel assignment. Based on the insights provided by our 
analytical model, we propose a simple decentralized algorithm 
which provides static channel assignments that are Nash equilibria 
in pure strategies for the objective of maximizing normalized 
network throughput. Our channel assignment algorithm requires 
neither any explicit knowledge of the topology nor any message 
passing, and provides assignments in only as many steps as there 
are channels. In contrast to prior work, our approach to channel 
assignment is based on the throughput metric. 

Index Terms — throughput modeling, fixed point analysis, chan- 
nel assignment algorithm, Nash equilibria, learning automata 



I. Introduction 




Fig. 1. A multi-cell infrastructure WLAN: DCF is used only for 
communications within the cells. Connectivity to nodes outside the 
cell is provided over a separate backbone which connects the APs to 
the Internet through a local server and a LAN router. 



This paper is concerned with infrastructure mode Wireless 
Local Area Networks (WLANs) that use the Distributed Co- 
ordination Function (DCF) Medium Access Control (MAC) 
protocol as defined in the IEEE 802.11 standard [1|. Such 
WLANs contain a number of Access Points (APs). Every 
client station (STA) in the WLAN associates with exactly one 
AP. Each AP, along with its associated STAs, defines a cell. 
Each cell operates on a specific channel. Cells that operate on 
the same channel are called co-channel. Thus, in our setting, 
DCF is used only for single-hop communication within the 
cells, and STAs can access the Internet only through their 
respective APs, which are connected to the Internet by a high- 
speed wire line local area network. Figure Q] depicts such a 
multi-cell infrastructure WLAN. 

To support the ever-increasing user population at high access 
speeds, WLANs are resorting to dense deployments of APs 
where, for every STA, there exists an AP close to the STA 
with which the STA can associate at a high Physical (PHY) 
rate [2|. However, as the density of APs increases, cell sizes 
become smaller and, since the number of non-overlapping 



channels is limitedj, co-channel cells become closer. Nodes 
in two closely located co-channel cells can suppress each 
other's transmissions via carrier sensing and interfere with 
each other's receptions causing packet losses. Thus, capacity 
might sometimes degrade with increased AP density 0, H. 
Clearly, effective planning and management are essential for 
achieving the benefits of dense deployments of APs. It has 
been demonstrated that a dense deployment of APs, along 
with careful channel assignment and user association control, 
can enhance the capacity by as much as 800% Q; but the 
technique has been tested on a small scald^- Large-scale 
WLANs are difficult to plan and manage since good network 
engineering models are lacking. 

Much of the earlier work on modeling WLANs deals with 
single- AP networks or the so-called single cells 0, 0, Q. 

'For example, the number of non-overlapping channels in 802. llb/g is 3 
and that in 802.11a is 12. 

2 In 0, to provide connectivity at high PHY rates, 24 APs were deployed 
in an area which could have been covered by a single AP at a low PHY rate. 



The existing performance analyses of multi-cell WLANs 
(i.e., WLANs consisting of multiple APs) are mostly based 
either on simulations or small-/medium-scale experiments 
El, E), 0. Studying large-scale WLANs by simulations 
and experiments is both expensive and time-consuming. It is, 
therefore, important to develop an analytical understanding 
of such WLANs in order to derive insights into the system 
dynamics. The insights thus obtained can be applied: (i) 
to facilitate planning and management, and (ii) to develop 
efficient adaptive schemes for making the WLANs self- 
organizing and self -managing. In this paper, we first develop 
an analytical model for 802.11 -based multi-cell WLANs and 
then apply our model to the task of channel assignment. 

Our Contributions: We make the following contributions: 
> We identify a condition, which we call the Pairwise 
Binary Dependence (PBD) condition, under which multi- 
cell WLANs can be modeled at the cell level (see Af3] in 
Section Hg. 

• We develop a scalable cell level model for multi-cell 
WLANs with arbitrary cell topologies (Section IVTi. 

« We extend the single cell TCP analysis of [8| to multiple 
interfering cells (Section [V-BI ). 

« We demonstrate how our model could be applied in 
conjunction with a Learning Automata (LA) algorithm 
for optimal channel assignment (Section [Villi) . 

• Based on the insights provided by our analytical model, 
we propose a simple decentralized algorithm which can 
provide fast static channel assignments (Section [DO . 

Unlike the node level models (9), iflOll or the link level 
models ifTTIl . lfl2l reported earlier in the context of multi- 
hop ad hoc networks, our cell level model does not require 
modeling the activities of every single node (resp. link). Thus, 
the complexity of our cell level model increases with the 
number of cells rather than the number of nodes (resp. links). 
However, our model can account for the number of nodes in 
each cell which may differ across the cells. 

We argue that the PBD condition is likely to hold in a 
dense deployments of APs (see Sections ITTIland irV- At . Hence, 
our cell level model can be applied to obtain a first-cut 
understanding of large-scale WLANs. 

Our cell level model is based on the channel contention 
model of Boorstyn et al. (9] and the transmission attempt 
model of [7|. Thus, our approach is similar to that of IflOll . 
However, our model of transmission attempt process is simpler 
(see the discussion following Equation O and the closed-form 
expressions for collision probabilities and cell throughputs that 
we derive are new. We also provide new insights. 

Our channel assignment algorithm provides assignments 
that are Nash equilibria in pure strategies for the normalized 
network throughput maximization objective (see Section llXb . 
Furthermore it provides an assignment in only M steps where 
M denotes the number of available channels. In contrast to 
prior work, our approach to channel assignment is based on the 
throughput metric (see the last paragraph in Section HJi. Our 



channel assignment algorithm is motivated by our analytical 
model which requires the knowledge of the topology of APs. 
However, our channel assignment algorithm requires neither 
any topology information nor any message passing. To be 
specific, our channel assignment algorithm takes the help of 
the network to discover an optimal assignment (in the sense 
mentioned above) in a completely decentralized manner. 

The remainder of this paper is organized as follows. In 
Section [TH we discuss the related literature. In Section Hill 
we provide the motivation for our simple cell level model. In 
Section llVl we provide our network model and summarize our 
key modeling assumptions. The analytical model is developed 
in Section [V] In Section I VII we validate our model by 
comparing with ns-2 simulations. In Section I VIII we apply 
our model to compare three different channel assignments for 
a 12-cell network and emphasize the impact of carrier sensing 
in 802.11 networks. In Section IVIIII we demonstrate how 
our analytical model could be applied along with a Learning 
Automata (LA) algorithm for optimal channel assignment. We 
propose a simple and fast decentralized channel assignment 
algorithm in Section IIX1 and summarize the conclusions in 
Section [X] Two lengthy derivations have been deferred to the 
appendices at the end of the paper. 

II. Related Literature 

Modeling of multi-hop ad hoc networks is closely related 
to that of multi-cell WLANs. In both cases, the hidden node 
and the exposed node problems lfl3ll are known to be the two 
main capacity-degrading factors. Multi-hop ad hoc networks 
and multi-cell WLANs are much harder to model than single 
cells precisely because of the presence of hidden and exposed 
nodes. Since the evolution of activities of each node is, in 
general, statistically different from that of the others', one 
needs to model the system at the node level J9), [ 10 1, or at the 
link level ifTTI . 1121 . i.e., the activities of every single node or 
link and the interactions among them need to be modeled. 

In the context of CSMA based multi-hop packet radio 
networks, Boorstyn et al. j£| proposed a Markovian model 
with Poisson packet arrivals and arbitrary packet length distri- 
butions. Wang and Kar [11] adopted the node level Boorstyn 
model to develop a link level model for 802.11 networks and 
studied the fairness issues. As a simplification, they assumed 
a fixed contention window. Garetto et al. flOl extended the 
Boorstyn model to multi-hop 802. 1 1 networks. They computed 
the packet loss probabilities and the throughputs per node 
accounting for the details of the 802.11 protocol. In particular, 
they incorporated the evolution of contention windows, which 
was missing in IfTTI . by applying the analysis of [7]. However, 
a direct application of the Boorstyn model required two 
different activation rates, namely, A and g, which complicated 
their model (see the discussion following Equation |2). We 
develop a simpler model by slightly modifying the Boorstyn 
contention model and making it particularly suitable for 802.1 1 
networks. 

In the context of multi-cell infrastructure WLANs, Nguyen 
et al. [ 14 1 proposed a model for dense 802.11 networks. To 



keep their interference analysis tractable, they assumed all 
the APs in the WLAN to be operating on the same channel. 
Recently, Bonald et al. ifTSIl applied the concept of "exclusion 
region" to model a multi-AP WLAN as a network of multi- 
class processor-sharing queues with state-dependent service 
rates. The concept of "exclusion domains" has also been 
applied in lfl2l to study the long-term fairness properties of 
large networks. The concept of exclusion says that, among a 
set of links that are in conflict, at most one can be active at 
any point of time. As argued in [12] and [15|, exclusion can 
be enforced in 802.11 networks by the RTS/CTS mechanism. 
However, suppression of interferers by RTS/CTS is not perfect 
because: (1) RTS frames can also collide and they must be 
retransmitted to enforce exclusion, (2) the overheads due to 
RTS/CTS that are transmitted at low rates and the capacity 
wastage due to RTS collisions might be significant especially 
when data payloads are transmitted at high rates, and most 
importantly, (3) RTS/CTS cannot completely eliminate hidden 
node collisions since nodes that cannot decode the CTS 
frames may collide with DATA frames that are longer than 
the Extended Inter Frame Space (EIFS). Thus, "exclusion" 
is strongly dependent on the virtual carrier sensing enforced 
by the RTS/CTS mechanism and is, essentially, a modeling 
approximation which ignores the possibility of hidden node 
collisions. The ineffectiveness of RTS/CTS is well known |[T6l 
and Basic Access can outperform RTS/CTS at high data rates 
even when hidden nodes are present ifPTl . 

Channel assignment in WLANs has been extensively stud- 
ied (see, e.g., fllfl, ED, (20|, ED, E), (23) and the 
references therein). Much of the existing work on channel 
assignment proposes to minimize the global interference power 
or maximize the global Signal to Noise and Interference Ratio 
(SINR) without taking into account the combined effect of the 
PHY and the MAC layers. Due to carrier sensing, nodes in 
802.1 1 networks get opportunity to transmit for only a fraction 
of time, and this must be accounted for when computing the 
global interference power or SINR. Such an approach is found 
only in |[T8"1 where the authors propose to maximize a quantity 
called "effective channel utilization". In reality, however, end 
users are more interested in the "throughputs" and, ideally, the 
objective should be to maximize the sum utility of throughputs. 
A heuristic method for computing throughputs can be found in 
[23 1 where the authors assume long term fairness among the 
nodes that are mutually interfering. The heuristic adopted in 
ll23l says that, ignoring collisions and channel errors, a node 
obtains approximately — of the maximum capacity where n 
denotes the cardinality of the maximum c%i«jl to which the 
node belongs. In reality, the size of the maximum clique to 
which a node belongs can be very different from the size of 
the maximum clique to which one of its interfering neighbor 
belongs, and this leads to inconsistent predictions for node 
throughputs. Nevertheless, the authors in ll23l claim that, their 
heuristic can predict the global throughput quite well even 
though individual node throughputs may not be. We remark 

3 A clique is a set of mutually interfering nodes. 



that, the multi-cell model in [15] does not actually model the 
DCF contention and adopts a model of capacity sharing similar 
to the heuristic of l23l . 

In summary, due to the inherent complexity of the problem, 
simplifying assumptions have often been made to develop 
tractable analytical model. The common simplification is to 
ignore hidden node collisions through suitable assumptions. 
We also ignore hidden node collisions to develop our model. 
However, unlike 11151 . our model can be applied either with 
the Basic Access mechanism or with the RTS/CTS mechanism 
provided that the PBD condition holds (see Af3]in Section llVl ). 
To our knowledge, an accurate throughput based approach 
which accounts for PHY-MAC interactions does not exist. 
Thus, we first develop a simple and accurate throughput model 
for multi-cell WLANs and then apply our model to the task 
of channel assignment. 

III. Motivation for a Cell Level Model 

In a dense deployment of APs with denser user population, 
it seems practically impossible to apply a node or a link level 
model for planning and managing the network. However, we 
can exploit a specific characteristics of dense deployments. 
Let R denote the cell radius, i.e., R is the maximum distance 
between an AP and the STAs associated to it. Let R cs denote 
the carrier sensing range, i.e., R cs is the distance up to which 
carrier sensing is effectiv^. We observe that, R cs >> R is 
likely to hold in a dense deployment of APs where, for every 
STA, there is an AP very close to the STA. With R cs » R, 
the network model can be simplified in the following ways: 

1) Since any transmitter 'T' is within a small distance 
R from its receiver 'R', a node 'H' that is beyond a 
distance R cs from 'T' (i.e., a hidden node) is unlikely 
to interfere with 'R^ i.e., carrier sensing would avoid 
much of the co-channel interference and we may ignore 
collisions due to hidden nodes. 

2) If Node- 1 in Cell-1 can sense the transmissions by Node- 
2 in Cell-2, then it is likely that all the nodes in Cell-1 
can sense the transmissions from all the nodes in Cell-2 
and vice versa, i.e., we may assume that nodes belonging 
to the same cell have an identical view of the rest of the 
network and interact with the rest of the network in an 
identical manner. Figure [2] summarizes this idea. 

The assumption that the AP and all its associated STAs have 
an identical view of the network has been applied in a dense 
AP setting [2| where the authors approximate STA statistics 
by statistics collected at the APs for efficiently managing their 
network. We adopt this idea of [2| to develop an analytical 
model. We identify the locations of the STAs with the locations 
of their respective APs, and treat a cell as a single entity, 
thus yielding a scalable cell level model. Unlike a node level 
model, the complexity of a cell level model increases with 
the number of cells rather than number of nodes. A simple 

4 Precise definition of carrier sensing range can be found in [13] and |24|. 
5 Ignoring noise, we have, SINR > , where v is called the path 

loss exponent and takes a value between 2 to 4. 




Fig. 2. A dense AP network: Small circular areas represent cells with 
radius R. APs are shown as tiny circles at the centers of the cells. The 
big circular disk represents the area covered by the carrier sensing 
range R CB of one of the APs. Observe that, with R ca >> R, other 
cells are either almost completely covered or almost not covered at all 
through carrier sensing. This indicates that, with R, cs >> R, nodes 
belonging to the same cell "perceive" identical medium activities, at 
least approximately. 

cell level model is particularly suitable for the task of channel 
assignment since channels are assigned to cells rather than to 
nodes. The foregoing simplification is extremely useful at the 
network planning stage when the locations of the STAs are not 
known but the locations of the APs and the expected number 
of users per cell might be known. Furthermore, since much of 
the traffic in today's WLANs is downlink, i.e., from the APs 
to the STAs, a large fraction of channel time is occupied by 
transmissions from the APs. It is then reasonable to develop a 
model based only on the topology of the APs and the expected 
number of users per cell assuming that the users are located 
close to their respective APs. 

The appropriate locations for placing the APs can be 
decided using standard RF tools and the physical topology 
of the APs can be determined. Given the physical topology 
of the APs and a channel assignment, we can obtain the 
logical topology, i.e., the topology of the co-channel APs 
corresponding to every channel. Given a logical topology, our 
analytical model can predict the cell throughputs accounting 
for the expected number of users per cell. The cell throughputs 
thus obtained can provide the goodness of the assignment to 
a channel assignment algorithm based on which the algorithm 
can determine a better channel assignment. Given the logical 
topology with the new channel assignment, our model can 
provide the modified feedback and this procedure can be 
repeated several times to arrive at an optimal assignment. We 
explain this approach in Section I Villi where we apply our 
model in conjunction with a Learning Automata (LA) based 
algorithm for optimal channel assignment. 

IV. Network Model and Assumptions 

We consider scenarios with R cs » R. Based on our 
earlier discussion, we assume that simultaneous transmissions 



by nodes that are farther than R cs from each other results 
in successful receptions at their respective receivers. We also 
assume that simultaneous transmissions by nodes that are 
within R cs always lead to packet losses at their respective 
receivers, i.e., we ignore the possibility of packet captur^. 
For simplicity, we assume symmetry among the nodes and 
of the propagation environment. In particular, this means that 
(i) nodes use identical protocol parameters (e.g., contention 
windows, retry limit etc.), and (ii) if Node-1 can sense Node- 
2, then Node-2 can also sense Node-1. 

We say that two nodes are dependent if they are within 
R cs ; otherwise, the two nodes are said to be independent. 
Two cells are said to be independent if every node in a cell 
is independent w.r.t. every node in the other cell; otherwise, 
the two cells are said to be dependent. Two dependent cells 
are said to be completely dependent if every node in a cell 
is dependent w.r.t. every node in the other cell. In this broad 
setting, our key assumptions are the following: 
A.l Only non-overlapping channels are used. 
A. 2 Associations of STAs with APs are static. This implies 

that the number of STAs in a cell is fixed. 
A. 3 Pairwise Binary Dependence (PBD): Any pair of cells 

is either independent or completely dependent. 
A.4 The STAs are so close to their respective APs that packet 

losses due to channel errors are negligible. 
A.5 The EIFS deferral has been disabled in the sense that 
medium access always starts after a Distributed Inter- 
Frame Space (DIFS) deferral. 

A. Discussion of the PBD Condition and Assumption A\5\ 

The PBD condition is crucial to the analytical model being 
developed in this paper. It includes the possibility that two 
given co-channel cells can be independent, i.e., two co-channel 
cells can be so far apart that activities in one cell do not 
affect the activities in the other cell. However, if two cells 
are dependent, the PBD condition rules out the possibility 
that only a subset of nodes in one cell can sense a subset of 
nodes in the other cell. We note that, in a dense deployment 
of APs, due to small cell radius R and R cs >> R, complete 
dependence would hold up to R cs and independence would 
hold beyond R cs , i.e., the PBD condition would hold, at least 
approximately. The PBD condition is a geometric property 
that enables modeling at the cell level since, the PBD 
condition holds, the relative locations of the nodes within a 
cell do not matter. Also, any interferer T, which must be 
within R cs of a receiver 'R', will also be within R cs of the 
transmitter T if the PBD condition holds. Hence, either (i) 
T can start transmitting at the same time as "F causing 
synchronous collisions at 'R', or (ii) T gets suppressed by 
T's RTS (or DATA) transmission followed by R's CTS (or 
ACK) transmission since CTS and ACK frames are given 
higher priority through SIFS « DIFS < EIFS). Thus, if 
the PBD condition holds, nodes do not require deferral by 

6 The difficulty in modeling packet capture lies in analytically computing 
the capture probabilities. Given the capture probabilities, packet capture can 
also be incorporated along the lines of 1251 . 



EIFS; deferral by DIFS would suffice. Thus, we assume that 
contention for medium access always begins after deferral by 
DIFS and we do not model the impact of EIFS. We remark 
that certain commercial products, e.g., Atheros cards, have 
user-configurable EIFS and EIFS can be easily disabled in 
such devices. 

V. Analysis of Multi-Cell WLANs with Arbitrary 
Cell Topology 

In this section, we develop a cell level model for WLANs 
that satisfy the PBD condition. We provide a generic model 
and demonstrate the accuracy of our model by comparing 
with simulations of specific cell topologies pertaining to: (i) 
networks with linear or hexagonal layout of cells (see Figures 
|3(a)|3(c)] i, and (ii) networks with arbitrary layout of cells 
(see Figure |3(d)| i. We index the cells by positive integers 
1, 2, . . . , TV, in some arbitrary fashion where N denotes the 
number of cells. Two completely dependent co-channel cells 
are said to be neighbours. Note that, two completely dependent 
cells operating on different non-overlapping channels are not 
neighbors. Let Af = {1,2, ... ,N} denote the set of cells 
and jVj (c Af) denote the set of neighboring cells of Cell-i 
(i S Af). Note that % £ M- Given the location of the APs and a 
channel assignment, we can obtain the set of neighboring cells 
for each cell. The key to modeling the cell level contention 
is the cell level contention graph Q which is obtained by 
representing every cell by a vertex and joining every pair 
of neighbors by an edge. Figures |3(a)|3(d)] also depict the 
contention graphs corresponding to each cell topology. 

In Section IV-AI we model the case where nodes are in- 
finitely backlogged and are transferring packets to one or more 
nodes in the same cell using UDP connection(s). Notice that 
single-hop direct communications among the STAs in the same 
cell without involving the AP are also allowed. In Section [V-BI 
we extend to the case when STAs download long files through 
their respective APs using persistent TCP connections. 

A. Modeling with Saturated MAC Queues 

Due to the PBD condition, nodes belonging to the same cell 
have an identical view of the rest of the network. When one 
node senses the medium idle (resp. busy) so do the other nodes 
in the same cell and we say that a cell is sensing the medium 
idle (resp. busy). Since the nodes are saturated, whenever a cell 
senses the medium idle, all the nodes in the cell decrement 
their back-off counters per idle back-off slot that elapses in 
their local medium^ and we say that the cell is in back-off. 
If the nodes were not saturated, a node with an empty MAC 
queue would not count down during the "medium idle" periods 
and the number of contending nodes would be time-varying. 
With saturated AP and STA queues, the number of contending 
nodes in each cell remains constant. 

We say that a cell transmits when one or more nodes in 
the cell transmit(s). When two or more nodes in the same 
cell transmit, an intra-cell collision occurs. Consider Figure 

7 Nodes belonging to different co-channel cells can have different views of 
the network activity. 
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Fig. 3. Examples of multi-cell systems: (a) four linearly placed co-channel 
cells, (b) five linearly placed co-channel cells, (c) seven hexagonally placed 
co-channel cells, and (d) seven co-channel cells with an arbitrary cell topology. 
The cell level contention graphs have also been shown where the dots represent 
the cells. Neighbors have been joined by edges. For example (a), the pairs 
{1,2}, {2,3} and {3,4} are dependent and the pairs {1,3}, {1,4} and 
{2, 4} are independent. 



|3(a)| There are periods during which all the four cells are in 
back-off. We model these periods, when none of the cells is 
transmitting, by the state $ where <£> denotes the empty set. The 
system remains in State <£> until one or more cell(s) transmit(s). 
When a cell transmits, its neighbors sense the transmission 
after a propagation delay and they defer medium access. We 
then say that the neighbors are blocked due to carrier sensing. 
However, two neighboring cells can start transmitting together 
before they could sense each other's transmissions resulting in 
synchronous inter-cell collisions. 

We observe that a cell can be in one of the three states: 
(i) transmitting, (ii) blocked, or (iii) in back-off. Modeling 
the synchronous inter-cell collisions requires a discrete time 
slotted model. However, this would require a large state space 
since the cells change their states in an asynchronous manner. 
For example, consider Figure |3(a)| and suppose that Cell- 
1 starts transmitting and blocks Cell-2 after a propagation 
delay. However, Cell-3 is independent of Cell-1 and can start 
transmitting at any instant during Cell-l's transmission. Thus, 
the evolution of the system is partly asynchronous and partly 
synchronous. To capture both, we follow a two-stage approach 
along the lines of 0~QTl . In the first stage, we ignore inter-cell 
collisions and assume that blocking due to carrier sensing is 
immediate. We develop a continuous time model as in |9j to 
obtain the fraction of time each cell is transmitting/blocked/in 
back-off. In the second stage, we obtain the fraction of slots 
in which various subsets of neighboring cells can start trans- 
mitting together. This would allow us to compute the collision 
probabilities accounting for synchronous inter-cell collisions. 
We combine the above through a fixed-point equation and 
compute the throughputs using the solution of the fixed-point 
equation. We define the following (7): 



n, := number of nodes in Cell-i 

/Si := (transmission) attempt probability (over the back-off 

slots) of the nodes in Cell-i 
ji := collision probability as seen by the nodes in Cell-i 

(conditioned on an attempt being made) 

Let Ai (T) denote the cumulative number of attempts made 
by a tagged node (any node due to symmetry) in Cell-i up to 
time T. Let Cj (T) denote the cumulative number of collisions 
as seen by the tagged node up to time T. Let B»(T) denote the 
cumulative number of back-off slots elapsed in Cell-i up to 
time T. Then, the attempt probability ft and the (conditional) 
collision probability 7$ of the tagged node are more precisely 
defined by 



ft 



lim 



A CO 



7i 



tto Bi(T) ' " : A, A,{T) ' 

Note that ft is the attempt probability of the nodes in Cell- 
i, irrespective of whether the nodes in the other cells can 
also attempt. This is a simplification and can be viewed as 
an extension to the decoupling approximation introduced in 
151 . Using the decoupling approximation and the analysis in 
1171 , the attempt probability ft of the nodes in Cell-i, Vi E Af, 
can be related to 7^ as 



Ci(T) 



ft = G(7<) 



1 + 7i . . . + 7. 



(1) 



60 + 7<6i • • • + 7*6* + • • • + 7f &x 

where K denotes the retry limit and bk, < k < K, denotes 
the mean back-off sampled after k collisions. 

The First Stage: When Cell-i and some (or all) of its neigh- 
boring cells are in back-off their (cell level) attempt processes 
compete until one of the cells, say, Cell-j, j £ Afi U {i}, trans- 
mits. Since, we ignore inter-cell collisions in the first stage, 
the possibility of two or more neighboring cells attempting 
together is ruled out. When Cell-i wins the contention, we say 
that it has become active. When Cell-i becomes active, it gains 
the control over its local medium by immediately blocking its 
neighboring cells that are not yet blocked. We assume that the 
time until Cell-i goes from the back-off state to the active state 
is exponentially distributed with mean -L. The activation rate 
Ai is given by 



A, 



1 - (1 - ft) r 



(2) 



where a denotes the duration of a back-off slot (in seconds) 
and 1 — (1 — j3i) ni is the probability that there is an attempt 
in Cell-i per back-off slot. Notice that we have converted 
the aggregate attempt probability in a cell per back-off slot 
to an attempt rate over back-off time. Also notice that, our 
assumption of exponential "time until transition from the back- 
off state to the active state" is the continuous time analogue of 
the assumption of geometric "number of slots until attempt" 
in the discrete time models of |6| and Q. 

Discussion 5.1: In [101, the authors use an unconditional 
activation rate A over all times as well as a conditional 



activation rate g over the back-off times and relate the two 
rates through a throughput equation which makes their model 
unnecessarily complicated. The activation rate A in our model 
is conditional on being in the back-off state. Thus, we use a 
single activation rate and our model is much simpler than that 
of (T0 |. Our modified approach can also be applied to simplify 
the node level model of flOl . ■ 
When Cell-i becomes active, it remains active and its neigh- 
bors remain blocked due to Cell-i until Cell-i's transmission 
finishes and an idle DIFS period elapses. The active periods of 
Cell-i are of mean duration — . When Cell-i becomes active 

Mi 

through a successful transmission (resp. an intra-cell collision) 
its neighbors remain blocked due to Cell-i for a success time 
T s (resp. a collision time T c fl Hence, — is given by 



1 

Ms 



m-l 



nMl - ft) ■ 
l-(l-ft)- )' {S 

1 - (1 - ft)™* 



{T c ) (3) 



where 



nift(l-ft) r 

1 - (i - ft) r 

irough i 



is the probability that Cell-i becomes 



active through a success given that it becomes active. 

Due to carrier sensing, at any point of time, only a set 
A (C Af) of mutually independent cells can be active together, 
i.e., A must be an independent set (of vertices) of the 
contention graph 60. From the contention graph Q, we can 
determine the set of cells Ba that get blocked due to A, and 
the set of cells Ua that remain in back-off, i.e., the set of 
cells in which nodes can continue to decrement their back-off 
counters. Note that A, Ba and Ua form a partition of Af, i.e., 
A, Ba and Ua are pairwise disjoint and Af = A U Ba U Ua- 
We take A(t), i.e., the set A of active cells at time t, as the 
state of the multi-cell system at time t. It is worthwhile now 
to mention the insensitivity result of Boorstyn et al. [9| which 
says that the product-form solution provided by their model 
is insensitive to the packet length distribution and depends 
only on the mean packet lengths. Applying their insensitivity 
argument, we take the active periods of Cell-i to be i.i.d. ex- 
ponential random variables with mean — . Due to exponential 
activation rates and active periods, at any time t, the rate of 
transition to the next state is completely determined by the 
current state A(t). For example, Cell-j, j S Ua, joins the set 
A (and its neighboring cells that are also in Ua join the set 
Ba) at a rate Xj. Similarly, Cell-i, i g A, leaves the set A (and 
its neighboring cells that are blocked only due to Cell-i leave 
the set Ba) to join the set Ua at a rate /i^. In summary, the 
process {A(t), t > 0} has the structure of a Continuous Time 
Markov Chain (CTMC). This CTMC contains a finite number 
of states and is irreducible. Hence, it is stationary and ergodic. 

8 For the Basic Access (resp. RTS/CTS) mechanism, T s corresponds to the 
time DATA-SIFS-ACK-DIFS (resp. RTS-SIFS-CTS-SIFS-DATA-SIFS-ACK- 
DIFS) and T c corresponds to the time DATA-DIFS (resp. RTS-DIFS). When 
DATA payload sizes are not fixed, T s and T c are to be computed using the 
expected payload size. 

9 An independent set of vertices of a graph G is a set of vertices of G such 
that no two vertices in the set are connected by an edge. 




Fig. 4. The CTMC describing the cell level contention for the four 
linearly placed cells given in Figure [3(a)] 



The set of all possible independent sets which constitutes the 
state space of the CTMC {A(t),t > 0} is denoted by A. 
For a given contention graph, A can be determined. For the 
topology given in Figure [3(a)] we have TV = {1,2,3,4} and 
A = {$,{1}, {2}, {3}, {4}, {1,3}, {1,4}, {2, 4}} where we 
recall that $ denotes the empty set. The CTMC {A{t),t > 0} 
corresponding to this example is given in Figure |4] 

It can be checked that the transition structure of the CTMC 
{A(t),t > 0} satisfies the Kolmogorov Criterion for re- 
versibility (see [26]). Hence, the stationary probability distri- 
bution ir (A), A £ A, satisfies the detailed balance equations, 

n(A)\i = tt(A U {i})p l , (Vi £ U A ) 
and the stationary probability distribution has the form 



■>ieA J 



(4) 



where pi :— ^ and 7r(<I>) (which denotes the stationary 
probability that none of the cells is active) is determined from 
the normalization equation 



E *w = L 

AeA 

Combining Equations [4] and [5] we obtain 



(5) 



n(A) 




En 



-, (VAeA) 



(6) 



Pj 



AeA ye- 4 



Convention: A product Y[ over an empty index set is taken to 
be equal to 1 . Notice that, with this convention, Equation [6] 
holds for 7r($) as well. 

The Second Stage: We now compute the collision proba- 
bilities 7j's accounting for inter-cell collisions. Note that ji 



is conditional on an attempt being made by a node in Cell- 
i. Hence, to compute 7,, we focus only on those states in 
which Cell-i can attempt. Clearly, Cell-i can attempt in State- 
A iff it is in back-off in State- A i.e., iff i £ Uj^. In all 
such states a node in Cell-i can incur intra-cell collisions 
due the other nodes in Cell-i. Furthermore, some (or all) of 
Cell-i's neighbors might also be in back-off in State- A If 
a neighboring cell, say, Cell-j, j £ TV, is also in back-off 
in State-A i.e., if j £ then a node in Cell-i can incur 
inter-cell collisions due to the nodes in Cell-j. The collision 
probability 7, is then given by Equation[7](appears at the top of 
the next page). A formal derivation of Equation [7] is provided 
in Appendix lAl 

Fixed Point Formulation: Equations[2][3]|6][7]and pi := — 
can express the 7/s as functions of only the p/s. Together 
with Equation Q] they yield an TY-dimensional fixed point 
equation where we recall that N is the total number of cells. 
The A^-dimensional fixed point equation can be numerically 
solved to obtain the collision probabilities 7j's, the attempt 
probabilities /3i's and the stationary probabilities 7r(„4)'s. In a ll 
of the cases that we have considered, the fixed point iterations 
were observed to converge to the same solutions irrespective 
of starting points. However, we have not yet been able to 
analytically prove the uniqueness of solutions. 

Calculating the Throughputs: The stationary probabilities 
of the CTMC {A(t),t > 0} can provide the fraction of time 
Xi for which Cell-i is unblocked. A cell is said to be unblocked 
when it belongs to either A or Ua- Thus, Vi G N, 



Xi= E n{A). 
AeA : ieAuUA 



(8) 



Definition 5.1: Let Qi, i e TV, denote the subgraph obtained 
by removing Cell-i and its neighboring cells in Mi from the 
contention graph Q. For a given contention graph Q, let A be 
defined as follows: 



a = e n 



(9) 



AeA ye- 4 



Let A; denote the A corresponding to the subgraph Qi. ■ 
An important observation which facilitates the computation 

of the Xi's is given by the following theorem. 

Theorem 5.1: The fraction of time Xi for which Cell-i is 

unblocked is given by 



(1 



A 



Vi £TV. 



(10) 



Proof: See Appendix iBl ■ 
Let 6; denote the aggregate throughput of Cell-i in a given 
multi-cell network and let & ni ,singieceii denote the aggregate 
throughput of Cell-i if it was an isolated cell containing 
iii nodes. Both 8^ and Q ni ,singieceU are in packets/sec. We 
approximate 0; by 



It = 



ieu A 



tt(A) 



(Vt e TV) 



(7) 



e, : = 



•^i ' ®ni,singlecell 
(l+/5 l )A 



A 



®rii,singlecell ? 



(ID 



in 



and 0j divided by rij gives the per node throughput 
Cell-i, i.e., $i = si (packets/sec). 

Discussion 5.2: Equation QT| is justified as follows. If Cell- 
i is indeed independent of every other cell in the network 
it is never blocked and does not incur inter-cell collisions. 
Then, we have Xi — 1 and O, = O 



rii ,singlecell 



However, in 



general, Cell-z gets blocked due to its neighbors for a fraction 
of time 1 — Xi and remains unblocked for a fraction of time 
Xi. If we ignore the time wasted in inter-cell collisions, the 
times during which Cell-z is unblocked would consist only 
of the back-off slots and the activities of Cell-i by itself. 
Thus, we approximate the aggregate throughput of Cell-i, over 
the times during which it is unblocked, by Q ni ,singleceU an d 
®ni,singlecell multiplied with x\ gives the aggregate through- 
put Qi of Cell-i in the multi-cell network. ■ 

Clearly, Equation [TT] is an approximation since the time 
wasted in inter-cell collisions have been ignored. However, 
we prefer to keep the approximation because: (1) it is quite 
accurate when compared with the simulations (see Section 
IVI1) . and (2) it can be efficiently computed since Q ni ,singlecett 
follows from a single cell analysis and the A as well as the 
Aj's can be computed using efficient algorithms [27 1. 

Complexity of the Model: In general, obtaining the state 
space A by searching for all possible independent sets A 
could be computationally expensive and the complexity grows 
exponentially with the number of vertices in the contention 
graph 11271 . For realistic topologies, where connectivity in the 
contention graph is related to distance in the physical network, 
efficient computation of A is possible up to several hundred 
vertices in the contention graph l27l . Thus, a cell level model 
is extremely helpful in analyzing large-scale WLANs with 
hundreds of cells since, unlike a node level model, each vertex 
in the contention graph now represents a cel0 

Large p Regime: Let rj (resp. rji) denote the number of 
Maximum Independent Seto (MISs) of Q (resp. Qi) (see 
Definition 15. U . Notice that, rji is also equal to the number 
of MISs of Q to which Cell-i belongs. From Equation [6] it is 
easy to see that, as pi — > oo, Vi G TV, we have, 



, ., - if A is an MIS, 

ir(A) < ?/ 

I otherwise. 

'"Note that, the state space A depends only on the cell topology and does 
not depend on the number of nodes in each cell. It also does not depend on 
whether the nodes use same/different protocol parameters. 

1 1 A maximum independent set of a graph is an independent set of the graph 
having maximum cardinality. 



Also, from Equations M IS] and [TJJ we observe that, as 

Pi — > oo, \/i G TV, we have, 

Hi 

V 

where we recall that Xi is the fraction of time for which Cell-i 
is unblocked. The quantity 



Qrii ,singlecell 



can be interpreted as the throughput of Cell-i, normalized 
with respect to Cell-i 's single cell throughput. We define the 

normalized network throughput Q by 



JV 



(12) 



Clearly, as pi — > oo, Vi G TV, the cells that belong to every 
MIS of the contention graph obtain normalized throughput 1 
and the cells that do not belong to any MIS obtain normalized 
throughput 0. Similar observations have also been made in 
[11]. We further observe that, as pi — > oo for all i G TV, only 
an MIS of cells can be active at any point of time. Since an 
MIS is always active, as pi — > oo, Vi G TV, we have, 

9 - a(Q) , 

where a(Q) denotes the cardinality of an MIS of Q and is also 
called the independence number of Q. Notice that a(Q) is a 
measure of spatial reuse in the network. 

B. Extension to TCP Traffic 

We now extend the analysis of Section IV-AI to the more 
realistic case when users access a local proxy server via 
persistent TCP connections. Our extension to TCP-controlled 
long file downloads is based on the single cell TCP-WLAN 
interaction model of [8|. The model proposed in [ 8 1 has been 
shown to be quite accurate when: (1) the local proxy server 
is connected with the AP by a relatively fast wired LAN such 
that the AP in the WLAN is the bottleneck, (2) every STA has 
a single persistent TCP connection, (3) there are no packet 
losses due to buffer overflow, (4) the TCP timeouts are set 
large enough to avoid timeout expirations due to Round Trip 
Time (RTT) fluctuations, and (5) the delayed ACK mechanism 
is disabled. We keep the above assumptions in this paper. 

In JSl, the authors propose to model a single cell having 
an AP and an arbitrary number of STAs with long-lived 
TCP connections by an "equivalent saturated network" which 
consists of a saturated AP and a single saturated STA. "This 
equivalent saturated model greatly simplifies the modeling 
problem since the TCP flow control mechanisms are now 
implicitly hidden and the total throughput can be computed 



using the saturation analysis |8)." Using the equivalent satu- 
rated model of [8|, the analysis of Section [V-AI can be applied 
to TCP-controlled long file downloads, taking = 2, Vi G Af. 
We explain the intuition behind this straightforward extension 
in the following. 

The equivalent saturated model of (8| is quite accurate in 
the single cell scenario because of the following reasons. For 
TCP-controlled long file transfers in WLANs, the AP has 
to serve n > 1 STA(s) by sending TCP DATA (resp. TCP 
ACK) packets to na downloading (resp. n u uploading) STAs 
(ni + n u = n). However, since DCF is packet fair, the AP 
does not get any prioritized access to the medium. Hence, 
in steady state, the AP's MAC queue never becomes empty, 
i.e., the AP is saturated"^!. Furthermore, since it is assumed 
that there are no buffer losses and no TCP timeouts, initially, 
the TCP windows keep growing and, in steady state (i) the 
TCP windows stay at their maximum, say, W max , for every 
connection, and (ii) the sum of the number of packets in the 
AP and in the STAs is equal to nW max . Most of these packets 
reside in the AP queue due to the DCF contention mentioned 
above and only a few STA queues remain non-empty at any 
point of time. It turns out that, the mean number of STAs 
with non-empty MAC queues w 1.5 for sufficiently large n 
El (n > 5 sufficesfl. Thus, the n STAs can be approximated 
by a single saturated STA and the equivalent saturated model 
can still predict the throughputs quite accurately. In fact, the 
equivalent saturated model can also accurately predict the 
collision probability of the AP which is indeed saturatec0 

Our extension to multiple cells is based on our key obser- 
vations that, 

1. As in single cells, the AP continues to be saturated in 
multi-cell scenarios as well. 

2. The number of non-empty STAs changes only at the end 
of successful transmissions 1281 : when the AP (resp. a 
STA) succeeds, this number increases (resp. decreases) 
by one. 

3. Whenever a cell is blocked, its activities are frozen, and 
hence, the number of non-empty STAs cannot change 
during the blocked periods. Thus, we can strip off the 
blocked periods and put the unblocked periods together 
over which time the cell would behave like a single cell. 

Applying the third observation above, the TCP download 
throughput of a cell, say, Cell-i, in the multi-cell scenario 
can be obtained by multiplying the single cell TCP download 
throughput of Cell-z with the fraction of time xi for which 
Cell-i is unblocked. Clearly, due to blocking, the likelihood 
of TCP timeouts increases. However, in our simulations, we 
have observed that, the minimum value of Retransmission 

12 This observation holds in steady state even for the n = 1 STA case if 
the maximum TCP window is large enough. In fact, for the n = 1 STA case, 
in steady state, the STA is also saturated |7 |. Thus, for the n = 1 STA case, 
the equivalent saturated model is exact. 

13 In fact, applying the Markov model of 1281 . it can be formally proved 
that the expected number of non-empty STAs is equal to 1.5. 

14 Our simulations indicate that, the equivalent saturated model cannot 
accurately predict the collision probabilities of the STAs. 



TABLE I 

Single Cell Results: Columns 2-3 Correspond to n Saturated 
Nodes. Columns 4-7 Correspond to TCP Downloads with n STAs 



n 


sat 

ry a a<. 

1 ana 


® ana 

(pkts/sec) 


' Sim 


AP 

i ana 




Sim 

(pkts/sec) 


ana 

(pkts/sec) 


1 





801.78 


0.0578 


0.0586 


454.01 


456.53 


2 


0.0586 


349.94 


0.0538 


0.0586 


456.06 


456.53 


3 


0.1077 


236.09 


0.0533 


0.0586 


456.09 


456.53 


4 


0.1473 


176.63 


0.0528 


0.0586 


456.17 


456.53 


5 


0.1812 


140.29 


0.0531 


0.0586 


456.05 


456.53 


6 


0.2100 


115.89 


0.0530 


0.0586 


456.10 


456.53 


7 


0.2348 


98.43 


0.0536 


0.0586 


456.88 


456.53 


8 


0.2565 


85.35 


0.0531 


0.0586 


456.97 


456.53 


10 


0.2927 


67.11 


0.0531 


0.0586 


456.02 


456.53 



TimeOut (RTO) is set to 200ms (which is the default value 
in ns-2.31), the analytical model is capable of predicting AP 
statistics quite well both in the single cell and the multi-cell 
scenarios. In fact, cells switch between the blocked and the 
unblocked states at a fast time scale of few packet transmission 
times which takes only a few milliseconds. Hence, timeout 
expirations are rare given that there are no buffer lossej^l 

VI. Results and Discussion 

We carried out simulations using ns-2.31 l29l . We created 
the example topologies given in Figure |3(a)|3(d)] We chose 
the cell radii and the distances among the cells such that the 
PBD condition holds. Nodes were randomly placed within the 
cells. The saturated case was simulated with high rate CBR 
over UDP connections. For the TCP case, we created one TCP 
download connection per STA. Each TCP connection was fed 
by an FTP source with the TCP source agent attached directly 
to the AP to emulate a local proxy server. The AP buffer 
was set large enough to avoid buffer losses. The EIFS deferral 
and the delayed ACK mechanism were disabled. Each case 
was simulated 20 times, each run for 200 sec of "simulated 
time". We report the results for "Basic Access". Similar results 
were obtained with "RTS/CTS". We took 1 1 Mbps data rate 
and packet payloads of 1000 bytes. The function 'fsolvef)" of 
MATLAB was used for solving the iV-dimensional fixed point 
equation. 

Table Q] summarizes the results for a single cell. Columns 
2-3 (resp. 4-7) correspond to the saturated case (resp. TCP 
download case) with n saturated nodes (resp. 1 AP and n 
STAs). The analytical results for the saturated case were 
obtained using [7] and that for the TCP case were obtained 
using (8). These single cell results obtained from known 
analytical models serve as the basis of our multi-cell results. 
The analytical throughputs per-node (resp. of AP) in Column 
3 (resp. Column 7) multiplied with the Xi's obtained from our 
multi-cell analysis provide the analytical throughputs per-node 
(resp. of AP) in the multi-cell cases (see Equation [TTJ. 

15 In certain cases, where some cells remain severely blocked, their blocked 
periods can be very long and TCP timeouts can occur for their respective 
connections. As we show in Section I VII our analytical model correctly 
identifies such severely blocked cells by predicting close to zero throughputs 
for them. 



TABLE II 

Results for t he Four Linearly Placed Cells Given in Figure 
|3(a)| when Each Cell Contains n = 5 Nodes 



Cell 






Q sim 


Sana 


Soo 


index 






(pkts/sec) 


(pkts/sec) 


(pkts/sec) 


1 


0.2351 


0.2399 


94.48 


97.41 


93.53 


2 


0.3005 


0.3146 


41.21 


46.66 


46.76 


3 


0.2999 


0.3146 


41.66 


46.66 


46.76 


4 


0.2359 


0.2399 


93.99 


97.41 


93.53 



TABLE III 

Results for t he Five Linearly Placed Cells Given in Figure 
|3(b)| when Each Cell Contains n = 5 Nodes 



Cell 
index 






S sim 

(pkts/sec) 


Sana 

(pkts/sec) 


Soo 

(pkts/sec) 


1 


0.1882 


0.1897 


129.35 


131.35 


140.29 


2 


0.3321 


0.3975 


8.69 


8.64 





3 


0.1892 


0.1925 


123.35 


126.41 


140.29 


4 


0.3321 


0.3975 


8.72 


8.64 





5 


0.1884 


0.1897 


129.31 


131.35 


140.29 



TABLE IV 

Results fo r the Seven Hexagonally Placed Cells Given in 
Figure [3(c)1 when Each Cell Contains n = 10 Nodes 



Cell 
index 


Tsim 


^ana 


S sim 

(pkts/sec) 


Sana 

(pkts/sec) 


Soo 

(pkts/sec) 


1 


0.2335 


0.8896 


0.003 


0.02 





2 


0.3045 


0.3158 


31.97 


32.35 


33.56 


3 


0.3061 


0.3158 


31.93 


32.35 


33.56 


4 


0.3055 


0.3158 


32.05 


32.35 


33.56 


5 


0.3054 


0.3158 


31.86 


32.35 


33.56 


6 


0.3070 


0.3158 


32.00 


32.35 


33.56 


7 


0.3058 


0.3158 


31.95 


32.35 


33.56 



TABLE V 

Resu lts fo r the Seven Arbitrarily Placed Cells Given in 
Figure [3(d)1 when CELL-i, 1 < i < 7 Contains n» = i + 1 Nodes 



Cell 


rii 






S sini 


Sana 


Soo 


index i 








(pkts/sec) 


(pkts/sec) 


(pkts/sec) 


1 


2 


0.0669 


0.0666 


320.66 


325.26 


349.94 


2 


3 


0.1163 


0.1163 


216.19 


219.65 


236.09 


3 


4 


0.2764 


0.3280 


12.48 


12.97 





4 


5 


0.3105 


0.3318 


34.29 


40.20 


46.76 


5 


6 


0.2505 


0.2585 


83.77 


84.92 


77.26 


6 


7 


0.3574 


0.3787 


28.67 


32.40 


32.81 


7 


8 


0.3062 


0.3139 


56.87 


59.21 


56.90 



A. Results for the Saturated Case 

Tables HIllIVI summarize the results for the example multi- 
cell cases depicted in Figures |3(a)|3(c)] respectively, when 
each cell contains n saturated nodes. Table FVl summarizes the 
results for the example case given in Figure [3(d)| when Cell-i, 
1 < i < 7, contains m = i + 1 saturated nodes. Quantities 
denoted with a subscript "sim" (resp. "ana") correspond 
to results obtained from ns-2 simulations (resp. fixed point 
analysis). In each case, 6^ represents the throughput per-node 
obtained by taking pi — > oo, Vi G M . We report only the mean 
values for our simulation results. The 99% confidence intervals 
were observed to be within 5% of the mean values. 

We show the plots corresponding to Table [V] in Figures 
[5] and [6] which compare the collision probability 7 and the 
throughput per node 6, respectively. In Figures and [6] we 
also show the relevant single cell results obtained from Table 
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Fig. 5. Comparing collision probability 7 for the example scenario 
in Figure [3(d)l when Cell-i, 1 < i < 7, contains n% = i + 1 saturated 
nodes. 
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Fig. 6. Comp aring throughput per node 8 for the example scenario 
in Figure [3(d)] when Cell-i, 1 < i < 7, contains m = i + 1 saturated 
nodes. 



H] i.e., the results one would expect had the seven cells been 
mutually independent. Referring Tables ITIirVl and Figures [5] 
and [6] we make the following observations: 
O-l.) Collision probabilities (resp. throughputs) in the multi- 
cell scenarios are always higher (resp. lower) than the corre- 
sponding single cell values (see Figures [5] and |6]l because (a) 
inter-cell collisions can be significant, and (b) due to inter-cell 
blocking, cells get opportunity to transmit only a fraction of 
time. 

0-2.) Our analytical model is quite accurate (less than 10% 
error in most cases) in predicting the collision probabilities 
and throughputs. However, our model always over-estimates 
the throughputs since the time wasted in inter-cell collisions 
have been ignored. Ignoring inter-cell collisions in the first 
stage of the model also over-estimates the fraction of time 
spent in back-off. Thus, the collision probabilities are also 
over-estimated. 

0-3.) The relative mismatch between the analytical model and 
the simulation is observed to be the worst for cells that remain 



blocked most of the time. For example, consider the first row 
of Table HVl which corresponds to Cell-1 in Figure [3(c)| Since 
Cell-1 is dependent with respect to every other cell, it obtains 
very few attempt opportunities in the simulations and its 
collision probability had to be averaged over very few samples. 
The same argument applies to the collision probability for 
Cell-3 in Figure El 

0-4.) Our analytical model correctly identifies the severely 
blocked cells, e.g., Cell-2 and Cell-4 in Table Hill Cell-1 
in Table IIVI and Cell-3 in Table [V] obtain significantly less 
throughputs than the other cells in the respective networks. 
Furthermore, our model works well with either equal or 
unequal number of nodes per cell. 

0-5.) Throughput distribution among the cells can be very un- 
fair even over long periods of time. Furthermore, introduction 
of a new co-channel cell can drastically alter the throughput 
distributions. For example, compare Tables HI] and Hill Cell-2 
and Cell-4 severely get blocked if Cell-5 is introduced to the 
four cell network given in Figure |3(a)| 

0-6.) The throughput of a cell cannot be accurately determined 
based only on the number of interfering cells. Consider, for ex- 
ample, Figure [3(d)| Cell-3 and Cell-4 each have two neighbors 
but their per node throughputs 8 are quite different (see Figure 
[6]). In particular, 64 > #3 even though = 4 < = 5. This 
is due to Cell-7 which blocks Cell-6 for certain fraction of time 
during which Cell-4 gets opportunity to transmit whereas Cell- 
1 and Cell-2 are almost never blocked and Cell-3 is almost 
always blocked due to Cell-1 and Cell-2. Thus, topology plays 
the key role and heuristic methods based only on the number 
of neighbors as in rt25l/ would fail to capture the individual 
cell throughputs. 

B. Results for the TCP Download Case 

Referring Columns 4-7 of Table Q] it can be seen that The AP 
statistics is largely insensitive to the number of STAs. Also, 
0-7.) The collision probability of the AP in a single cell with 
any number of STAs is approximately equal to the collision 
probability in a single cell with two saturated nodes. This 
can be verified by comparing Columns 4 and 5 with Row 2 
Column 2 in Table U 

0-8.) The AP throughput does not change with the num- 
ber of STAs. In fact, the AP throughput with any number 
of downloading STAs is equal to the per node throughput 
in a single cell containing two saturated nodes with pay- 
load size l tcp-data+Ltc P -ack where Ltcp-data and 
Ltcp~ack denote the size of TCP DATA and TCP ACK 
packets. 

Observations 0-7 and 0-8 are well-known and they form 
the basis for the equivalent saturated model of [8 1. Observation 
0-7 has led to the conclusion in 1 3 1 that the collision probabil- 
ity of the nodes in a WLAN containing m mutually interfering 
(i.e., dependent) APs is equal to the collision probability in 
a single cell containing 2m saturated nodes since each cell 
can be assumed to consist of two saturated nodes (a saturated 
AP and a saturated STA). Tables IVIHVIIII summarize the AP 
statistics for the topologies in Figures |3(a)| |3(b)| and |3(d)| 



TABLE VI 

Results for the AP Corresponding to Figure [3(a)| when Each 
Cell Contains 1 AP and n = 5 STAs 



Cell 


lsim,AP 


Tana, AP 


Vsim,AP 


&ana,AP 


0oo, AP 


index 






(pkts/sec) 


(pkts/sec) 


(pkts/sec) 


1 


0.1038 


0.1033 


306.33 


318.73 


304.35 


2 


0.1560 


0.1574 


153.16 


169.18 


152.18 


3 


0.1555 


0.1574 


153.06 


169.18 


152.18 


4 


0.1038 


0.1033 


306.41 


318.73 


304.35 



TABLE VII 

Results for the AP Corresponding to Figure [3(b)] when Each 
Cell Contains 1 AP and n = 5 STAs 



Cell 
index 


1sim,AP 


7ana,AP 


6sim,AP 

(pkts/sec) 


9ana,AP 

(pkts/sec) 


0oo, AP 

(pkts/sec) 


1 


0.0728 


0.0775 


381.21 


387.16 


456.53 


2 


0.1793 


0.1950 


75.16 


85.62 





3 


0.0744 


0.0832 


340.24 


346.47 


456.53 


4 


0.1786 


0.1950 


75.23 


85.62 





5 


0.0728 


0.0775 


381.15 


387.16 


456.53 



TABLE VIII 

Results for the AP Corresponding to Figure [3Cd)] when Each 
Cell Contains 1 AP and n = 5 STAs 



Cell 


~isim,AP 


7ana,AP 


Vsim,AP 


Vana,AP 


#00, AP 


index 






(pkts/sec) 


(pkts/sec) 


(pkts/sec) 


1 


0.0610 


0.0670 


421.70 


425.83 


456.53 


2 


0.0604 


0.0670 


421.92 


425.83 


456.53 


3 


0.2010 


0.2528 


33.79 


38.50 





4 


0.1561 


0.1685 


141.80 


156.41 


152.18 


5 


0.0987 


0.1028 


317.39 


329.06 


304.35 


6 


0.1551 


0.1644 


158.55 


172.64 


152.18 


7 


0.1061 


0.1099 


301.15 


314.10 


304.35 



respectively. Since the AP statistics does not change with the 
number of STAs, we report the results with n = 5 STAs in 
each case. We show the plots corresponding to Table lyTTTl in 
Figures [7j and [8] which compare the collision probability 7 
and the throughput per node 9, respectively. Referring Tables 
IVIHVIIII and Figures [7j and [8] we conclude that the foregoing 
observations (O.I-O.6) for the saturated case carry over to 
TCP-controlled long file transfers as well. Furthermore, we 
generalize the conclusion of [3| as follows: 
0-9.) The collision probability in a WLAN depends, not only 
on the number of interfering APs but also on the fraction 
of time for which the neighboring APs can cause collisions. 
In general, collision probability does not grow as twice the 
number of interfering APs and depends on the cell topology. 

Referring Tables IVIHVIIII we extend the validity of the 
equivalent saturated model of [8| as follows: 
O-10.) The equivalent saturated model of [8| proposed in the 
context of a single cell, preserves its desirable properties, i.e., 
it predicts the AP statistics quite well when extended to a 
multi-cell WLAN that satisfies the PBD condition. 

C. Variation with p 

Till now, we have discussed the results corresponding to 
the payload size of 1000 bytes. We now examine how the 
results vary with the payload size. Figures l9l and [Tol depict the 
variation with payload size of analytically computed collision 
probabilities and normalized cell throughputs, respectively, for 
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Fig. 7. Comparing collision probability 7 for the example scenario 
in Figure [3(d)] when each cell contains an AP and n = 5 STAs. STAs 
are downloading long files through their respective APs using TCP 
connections. 
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Fig. 9. Variation of collision probabilities with payload size for the 
seven cell network in Figure |3(d)| when each cell contains n = 10 
saturated nodes. 
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Fig. 8. Comparing throughput per node 6 for the example scenario 
in Figure [3(d)] when each cell contains an AP and n = 5 STAs. STAs 
are downloading long files through their respective APs using TCP 
connections. 



Fig. 10. Variation of normalized cell throu ghputs with payload size 
for the seven cell network in Figure |3(d)| when each cell contains 
n = 10 saturated nodes. 



the seven cell network of Figure [3(d)| when each cell contains 
n = 10 saturated nodes. Similar results were obtained for the 
TCP case as well. From Figures [9] and [lOl we observe that: 
O-ll.) The results are largely insensitive to the variation in 
payload size. Moreover, as the payload size increases, the 
normalized cell throughputs become closer to the normalized 
throughputs under the infinite p approximation, i.e., they 
become closer to x\ = x-2 = 1,%3 = 0,2:4 = x§ = 
\iX§ = £7 = |. Hence, except for very small payload 
sizes, the infinite p approximation can be expected to provide 
fairly accurate predictions. Furthermore, for sufficiently large 
payload sizes, we have, 9 w a(Q). 

VII. A Simple Design Example 

We consider a 12-cell network as shown in Figure QT] and 
apply our throughput model to compare three different channel 
assignments as shown in Figures |12(a)|[l"2(c)| Notice that the 
"logical" contention graphs for the three assignments in Fig- 




Fig. 11. A 12-cell network: Cell indices have been marked. Also 
shown is the cell level physical contention graph assumed for the 
network. 

ures |12(a)|[T2"(c)| are different from the "physical" contention 
graph of the network shown in Figure QT] which is based 
on the physical separation among the cells. Two physically 
dependent cells become logically independent if they are 
assigned different channels and the corresponding edge will 
be missing in the logical contention graph. Our objective is to 




determine the best, among the three assignments, in terms of 
the normalized network throughput (see Equation [T2l i and 
the fairness index reported in [ 30 1 . Considering fairness at the 
cell level, the fairness index J is given by [30 1 

*(E£i»?)' 

where we recall that Xi is the normalized throughput of Cell- 
i. Note that, < J < 1, and higher values of J imply more 
fairness. 

Observe that co-channel cells in the first assignment have 
topology as in Figure |3(a)| There are three subsystem of 
co-channel cells each consisting of four cells. The infinite 
p approximation predicts that, the two cells at the middle 
(resp. at the end) for each subsystem of co-channel cells 
would obtain normalized throughput « -| (resp. « |). The 
normalized network throughput for the first assignment is 6 
and the fairness index is 0.45. In the second assignment, 
each cell obtains normalized throughput « |. The normalized 
network throughput for the second assignment is 4 and the 
fairness index is 1. Thus, the fairness has attained its maximum 
in the second assignment. However, this assignment is not 
desirable since for this assignment is only | of in the 
first assignment. In the third assignment, each cell obtains nor- 
malized throughput « i. The normalized network throughput 
for the third assignment is 6 and the fairness index is equal 
to 1. Clearly, the third assignment is the best one among the 
three. 

Note that, an SINR based model is likely to prefer the first 
assignment to the third. However, interference by nodes within 
R cs can be taken care by carrier sensing, and hence, such 
interference should not be given undue weights. If required, 
cells that are very close to each other may be assigned the 
same channel to form "big" cells. Unlike a cell which covers 
a large area by a single AP and provides connectivity at low 
PHY rates to STAs that are located at the edge of the cell, 
a cluster of mutually interfering small co-channel cells each 
containing an AP to ensure connectivity at high PHY rates 
can provide high capacity and fairness. This emphasizes the 
impact of carrier sensing in 802.11 networks. 



VIII. Applying the Analytical Model along with a 
Learning Algorithm 

In the 12-cell network of Figure [TT] with 3 non-overlapping 
channels, we need to examine 3 12 possibilities to determine 
the optimal assignment. It is also not clear if the normalized 
network throughput for the 12-cell network can be actually 
more than 6 with some other assignment. Under the infinite 
p approximation, the maximum value of is equal to the 
maximum independence number over all possible logical con- 
tention graphs. Clearly, the problem of finding an assignment 
that maximizes is NP-hard and, in general, it is not desirable 
to examine all the possibilities. 

We now demonstrate how our analytical model could be 
applied along with a Learning Automata (LA) algorithm 
called the Linear Reward-Inaction (Lr-j) algorithm lf3~TTl to 
obtain optimal channel assignments. For every cell i S Af, 
we maintain an A/-dimensional probability vector ip i where 
M denotes the number of available channels. The 
algorithm proceeds in steps. Let k = 0,1, 2,... be the step 
indices. At each step k, a channel is selected for Cell-i, i 6 Af, 
according to the probability distribution 

where 

Pi,j(k) :— the probability that Channel- j is selected for 
Cell-z at step k. 

Let c(k) = (ci(k), C2(fc), . . . , Cjv(fc)^ be the channel 
assignment at step k, where 

Ci(k) := the channel selected for Cell-i at step k. 

Given a channel assignment c, we can obtain the logical 
contention graph Q(c) and applying our throughput model, 
we can compute the normalized throughput Xi(c) of every 
cell i with the assignment c. We define the sum utility for an 
assignment c by 

N 

U(c) :=^u(xi{c)), (13) 
i=i 



where u(-) is some suitably defined increasing concave func- 
tion. For a given u(-), we compute the sum utility at step k 
by 

N 

u(c(k)) = j2<xiwm, 

i=l 

and apply the Lr_j algorithm which consists of the following 
steps EH: 

1. Begin with p itj G (0, 1), Vi G Af, j = 1, 2, . . . , M, such 
that Pj ■ 1 = 1, where 1 is the M-dimensional vector 
with all components equal to 1. 

2. Update the probability vectors of Cell-i, Vi G Af, as 

Pl {k + 1) = Pl (k) + b U(c(k))(8 Ci{k) - Pi (k)) , 

where Sj denotes the M-dimensional probability vector 
with unit mass on Channel- j and < b < 1. 
The parameter b is called the learning parameter or the step- 
size parameter which determines the convergence properties 
of the algorithm; with smaller b convergence is slower but the 
algorithm may not converge to the desired optimum if b is not 
small enough [31 1. Notice that, to ensure non-negativity of the 
probability vectors after every update with any b G (0, 1), the 
sum utility U must satisfy < U < 1. To maximize 6 = 
Sili x i> we ta ^ e trie avera g e normalized network throughput 

1 N 

i=l 

as the sum utility U where we compute the X{'s applying the 
infinite p approximation. Note that, < Uq < 1. Other utility 
functions can also be considered provided that < U < 1. 
Define 

p : = (Pi,P 2 >---, Pa- 
using the terminology of BP . the super vector P will 
be called a strategy for the channel assignment problem. 
A strategy P such that, Vi G Af, pij — 1 for some j, 
j = 1,2, ... , M, will be called a pure strategy. A strategy P 
that is not a pure strategy is called a mixed strategy. Under the 
Lr-i algorithm, {P(fc), k > 0} is a Markov process with the 
pure strategies as the only absorbing states. Thus, in practice, 
the Lr_j algorithm always converges to a pure strategy rather 
than to a mixed strategy j3Tj. Furthermore, by Theorem 4.1 
of (|3Tj, the Lr_j algorithm always converges to one of the 
Nash equilibria. Thus, in practice, channel assignment by the 
Lit-i algorithm always provides an assignment c* which is 
one of the Nash equilibria in pure strategies in the sense that 

U(c*) > U(c), 

for all c that differs from c* by exactly one element, i.e., 
changing the channel of one of the cells in the assignment c* 
does not increase the sum utility U. 
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Fig. 13. Evolution of the average normalized network throughput 
Uq under the Lr—i algorithm for the 12-cell network in Figure [T7I 
with M = 3 channels, b = 0.01, and pij = i, Vi £ N and j = 
1,2,. ..,M. 

A. Results for Channel Assignment by the Lr-i Algorithm 

Figure [13] shows the evolution of the average normal- 
ized network throughput U® under the L#_j algorithm with 
M = 3 channels and b = 0.01 for the 12-cell network 
in Figure [TT] We began with uniform unbiased probability 
vectors p L = (i |, i), Vi G Af. The algorithm converged to 
the assignment jj c x = c 6 = c 9 = 1, c 4 = c 7 = c 12 = 3, c 2 = 
c 3 = C5 = c 8 = C10 = c n —2. Notice that Uq converges 
to I which corresponds to 9 = 8. It can be verified that 
this is the maximum possible value that 9 can take for the 
given network with 3 channelQ We found that the Lr^i 
algorithm always converges to a globally optimum solution if 
started with uniform unbiased probability vectors and if b is 
sufficiently small. However, starting with probability vectors 
that are highly biased towards some assignment or if b is 
not small enough, it may not converge to a globally optimum 
solution. We demonstrate this through Figures [T4lfT7l for the 
7 -cell example in Figure [3(d)| with M = 2 channels. 

Beginning with unbiased probability vectors and b = 0.01 
the algorithm may converge to a global optimum with U@ = 1 
or 9 — 7 (Figure [Pfl). There are two possible global optima, 
namely, ci = c% = C4 = C7 = 1, C3 = C5 = cq = 2 and 
ci = c 2 = c 4 = c 7 = 2, c 3 = c 5 = c 6 = 1. In Figure Q3] we 
show that, with b = 0.01, the algorithm may also converge 
to a solution that is not a global optimum. In fact, for the 
specific simulation run, the algorithm converged to the Nash 
equilibrium c\ = c 2 = C5 = cq = 1 , C3 = C4 = C7 = 2 
with 9 = 6 which is not a global optimum. The assignment 
ci = C2 = C5 = cq = 1, C3 = C4 = C7 = 2 is a Nash 
equilibrium in the sense that changing the channel of one of the 
cells does not increase 9, and hence, the sum utility Uq. With 

16 We do not provide the plots showing evolution of the Pi.j's to save space. 

17 Note that permutations of channels do not change the logical topology 
and the associated 0's would be equal. For example, ci = eg = eg = 
2, C4 = C7 = C12 = 1, C2 = C3 = C5 = eg = cio = cn = 3 would also 
give G = 8. 
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Fig. 14. Evolution of the average normalized network throughput C7@ 
for the 7-cell example in Figure [3(d)] with M = 2 channels, b = 0.01, 
and pn = pi2 = 0.5, % = 1, 2, . . . , 7. In this example, the Lb_-i 
algorithm converges to a global optimum. 
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Fig. 15. Evolution of the average normalized network throughput Uq 
for the 7-cell example in Figure [3(d)] with M = 2 channels, b = 0.01, 
and pn = pi2 = 0.5, i = 1, 2, . . . , 7. In this example, the Lr—i 
algorithm converges to a local optimum. 



Fig. 16. Evolution of the average normalized network throughput 
Uq for the 7-cell example in Figure [3(d)] with M = 2 channels, b = 
0.001, and pn = P21 = P32 = P42 = PBl = P61 = P72 = 0.8. In 
this example, the Lr_j algorithm converges to a global optimum. 
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Fig. 17. Evolution of the average norm alized network throughput 
Uq for the 7-cell example in Figure [3(d)] with M = 2 channels, b = 
0.001, and pn = p 2 i = P32 = P42 = Psi = Pel = P72 = 0.9. In 
this example, the Lr—i algorithm converges to a local optimum. 



b = 0.001, and beginning with unbiased uniform probability 
vectors, we observed that the algorithm always converged to a 
global optimum. With b = 0.001, the algorithm converged to a 
global optimum even when we began with probability vectors 
biased towards the assignment Ci = C2 = C5 = eg = 1, C3 = 
C4 = C7 = 2 (Figure [T6b. However, when the initial probability 
vectors were highly biased towards the assignment c\ = C2 = 
C5 = cq = 1 , C3 = C4 = cj = 2, the algorithm converged to 
the biased assignment which is not a global optimum (Figure 

us. 

IX. A Simple and Fast Algorithm for Maximizing 
Normalized Network Throughput 

We observe that the Lr—i algorithm takes a large number 
of iterations to converge and guarantees convergence only to 
Nash equilibria. Clearly, for the 7-cell example with M = 2 
channels, an exhaustive search would have required exam- 
ination of only 2 7 = 128 possibilities whereas the -Lr-j 



algorithm takes many more steps to converge. The Linear 
Reward-Penalty (Lr-p) algorithm of [21 1 and the simulated 
annealing algorithm of \22\ guarantee convergence to a glob- 
ally optimum solution as the number of iterations goes to 
infinity. A greedy version of simulated annealing algorithm 
in IF22I is relatively faster but still takes a large number 
of iterations to converge and guarantees convergence only 
to locally optimum solutions. To maximize the normalized 
network throughput 9 we now propose a simple and fast 
decentralized algorithm which can be easily implemented in 
real networks. 

We form a physical contention graph Q in which every 
completely dependent pair of cells are neighbors (since 
we have not yet assigned channels) and our objective is 
to transform Q to a logical contention graph Q{c) by an 
assignment c so that a(Q(c)) is maximized. As noted in 
Observation O.ll, for sufficiently large packet sizes, we have 
w a(Q). Hence, maximizing a(Q(c)) would maximize 0. 



Maximal Independent Set Algorithm (mISA) : We pro- 
pose the following channel assignment algorithm: (recall that 
M is the number of available channels) 

(1) Choose a maximal independent set of cells, 
assign them Channel- 1 and remove them from the graph. 

(2) Increment the channel index and repeat the procedure 
on the residual graph until one channel is left. 

(3) Assign Channel-M to all the cells in the residual graph 
after M — 1 steps. 

Notice that mISA takes only M steps. Also notice that, the 
residual graph may become null (i.e., it might not have any 
vertices left) after M' < M steps. Clearly, mISA is based 
on a classical graph coloring technique, but the novelty lies in 
our recognition of the notion of optimality that mISA provides 
which is given by the following 

Theorem 9.1: The channel assignments by mISA are Nash 
equilibria in pure strategies for the objective of maximizing 
normalized network throughput as pi — > oo, Vi £ Af. 
Furthermore, mISA provides a globally optimum solution if 
M > D(Q) + 1 where D{Q) denotes the maximum vertex 
degree in the physical contention graph Q. 

Proof: If the residual graph becomes null in less than M 
steps, then every cell would have a channel different from all 
its neighbors' and, we have, = N, i.e., mISA would provide 
a globally optimum solution. Hence, to prove the first part of 
the theorem, we assume, without loss of generality, that the 
residual graph after M — 1 steps is not null. 

Suppose that Nj cells are assigned Channel- j in Step-j, 1 < 
j < M — 1. Let Qj be the residual graph after j steps, 1 < j < 
M—l. Let Ok denote the aggregate normalized throughput of 
the cells on Channel-/c, 1 < fc < M. Then, we have, Qj = Nj, 
1 < j < M — 1, due to independence, and 6m = a(Gni-i)- 
Hence, = Ylk=i ®* = Y^jL^ 1 N j + «(6m-i)- Suppose 
that a cell on Channel-j, j ^ M, is moved to Channel-fc, 
k ^ j. Then, Qj decreases by 1 but Ok can increase by at 
most 1 . Hence, cannot increase. Suppose now that a cell on 
Channel- M is moved to Channel-j, 1 < j < M — 1. Clearly, 
any cell on Channel-Af is dependent (in the original graph Q) 
w.r.t. at least one of the Nj cells on Channel-j since the Nj 
cells that are already on Channel-j, 1 < j < M — 1, form an 
mIS. Hence, Qj does not change but Qm can only decrease. 
Hence, cannot increase by changing the channel of only one 
of the cells and the first part of the theorem is proved. 

For the second part of the theorem, notice that, every cell 
in the residual graph Qj after j steps, 1 < j < M — 1, must 
be dependent (in the original graph Q) w.r.t. at least one cell 
on each channel I, 1 < I < j — 1. This follows from the 
maximality of the independent sets of cells chosen in each of 
the first M — 1 steps. In particular, every cell in the residual 
graph Gm-i is dependent (in the original graph Q) w.r.t. at 
least one cell on each channel j, 1 < j < M — 1. Since 
different channels are assigned in each step, the vertex degree 

18 The lower case 'm' corresponds to "maximal" as opposed to the upper 
case 'M' which corresponds to "maximum". 



of every cell in Qm-i restricted to the cells in Gm-i alone 
would be < D(G) — (M — 1). If M > D{Q) + 1, then the 
cells in Gm-i would form an independent set and every cell 
would have a channel different from its neighbors'. Thus, if 
M > D(Q) + 1, then mISA would provide an assignment with 
= N, i.e., a globally optimum solution. ■ 
Implementation of mISA: mISA can be implemented in 
a decentralized manner as follows. APs sample random back- 
offs using a contention window W and contend for accessing 
the medium using Channel- 1. When an AP wins the contention 
it keeps transmitting broadcast packets separated by Short Inter 
Frame Space (SIFS) for some duration T » aW where we 
recall that a is the duration of a back-off slot. This emulates 
the infinite p situation since and AP after wining the contention 
does not relinquishes the control over its local medium. We 
had observed that, as pi — > oo, Vi £ Af, only the cells that 
belong to an MIS obtain non-zero normalized throughputs. 
But this holds only in an ensemble average sense. If an AP, 
after wining the contention, does not relinquishes the control 
over its local medium, in a particular sample path, an mIS 
of APs (which may not be an MIS) would grab the channel 
during 71 19 l. This is not surprising since with infinite p^'s, the 
CTMC {A(t),t > 0} becomes absorbing with the maximal 
independent sets of cells as the only absorbing states and we 
cannot expect the time average to be equal to the ensemble 
average. 

At time T, APs that could transmit consecutive broadcast 
packets stop contending until time (M — 1) X T and APs 
that remain blocked switch to Channel-2, sample fresh back- 
offs and keep contending until 2T and so on. APs that remain 
blocked throughout the duration (M—l) x T stick to Channel- 
M. Thus, in every time duration T, a maximal independent 
set of APs would be assigned a channel. Normal network 
operation can begin after time (M — 1) x T. Notice that, 
mISA does not require any knowledge of AP topology and runs 
in a completely decentralized manner without any message 
passing. In addition, if there is a central controller to which 
the APs can communicate, mISA can be repeated several times 
before normal network operation could begin. The central 
controller, which obtains the global view of the channel 
assignments, can choose the best among the solutions provided 
by mISA. In absence of centralized control, mISA can be 
invoked periodically. Thus, mISA can be easily implemented 
in real networks in a completely decentralized manner if the 
number of channels for every AP is the same and known. 
However, mISA requires loose synchronization among the APs 
which necessitates some message passing before mISA could 
begin. 

X. Conclusions and Future Work 

In this paper, we identified a Pairwise Binary Dependence 
(PBD) condition that allows a scalable cell level modeling of 
WLANs. The PBD condition is likely to hold at higher PHY 

"if W is large enough, the possibility of two dependent APs transmitting 
together can be ruled out. 



rates and denser AP deployments. We developed a cell level 
model both under saturation condition and for TCP-controlled 
long file downloads. Our analytical model was shown to 
be quite accurate, insightful and capable of comparing few 
design alternatives. Thus, we believe that our modeling frame- 
work is a significant step toward gaining "first-cut" analytical 
understanding of WLANs with dense deployments of APs. 
We demonstrated how our analytical model could be applied 
along with the Linear Reward-Inaction learning algorithm for 
optimal channel assignment. Based on the insights provided by 
our model, we also proposed a simple decentralized algorithm 
called mISA which can provide channel assignments that are 
Nash equilibria in pure strategies in only as many steps as 
there are channels. Although mISA is based on a standard 
graph coloring technique, its applicability to the considered 
problem is strengthen by and its practical implementability 
is guided by a rigorous analytical model developed in this 
paper. Developing simple and practical algorithms for general 
objective functions, other than maximizing the normalized 
network throughput, is a topic of our ongoing research. 
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Consider a tagged node in Cell-i. Let 



Appendix A 
Derivation of Equation [7] 



Ai (T) : number of attempts made by the tagged node up to time T 
Ci (T) : number of collisions as seen by the tagged node up to time T 
Af-(T) : number of attempts made by the tagged node in State-^4 up to time T 
Cf-(T) : number of collisions as seen by the tagged node in State-„4 up to time T 
(T) : number of back-off slots elapsed in Cell-i in State-,4 up to time T 



Then, we have 
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where the last but one step follows from the facts that: (i) the time spent in State-„4 up to time T is equal to ir(A)T and 
if Cell-i is in back-off in State- A then B^(T) = 7r ^ T , (ii) the second limit within the brackets in the numerator is the 
long-run fraction of back-off slots in State-„4 in which the tagged node attempts which is equal to irrespective of State-„4 
given that Cell-i is in back-off in State-A and (iii) the third limit within brackets in the numerator is the long-run fraction 
of attempts made by the tagged node in State-„4 that result in collisions and depends on the set of neighbors that are also in 
back-off in State-A 



Appendix B 
Proof of Theorem 15. II 

Solving the balance equations for the CTMC {A(t), t > 0} (Equation|4]l together with the normalization equation, we obtain 
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Equation [8] can now be expanded as 
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(since .4 n = $) 
(by Equation @J 
(by Equation [TO]) 
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where the fourth step follows from the facts that: (i) if Cell-i is active, then none of its neighbors can be active, and (i) if 
Cell-i is in back-off, then neither Cell-i nor any of its neighbors can be active. Observe now that, the independent sets of Qi 
are also independent sets of Q. Thus, A, can be obtained by restricting the summation in Equation|9]to only those independent 
sets of Q that do not contain Cell-i or any cell in A/i, i.e., 

AeA : An({i}uAfi)=<S> ^ ieA ' AeA : An{i}=*, AfW t =9 / 

where the second equality follows from the fact that union of two sets will be equal to the empty set iff each of the sets is an 
empty set. Applying Equation [18] in the last step of Equation [17] we obtain Equation [TO] 



